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ABSTRACT 

We re-estimate the peculiar velocity of the Sun with respect to the local standard of rest 
(LSR) using a sample of local stars within 600 pc of the Sun, selected from the Large Sky 
Area Multi-Object Fiber Spectroscopic Telescope (LAMOST; also named the Guoshoujing 
Telescope) Spectroscopic Survey of the Galactic Anti-centre (LSS-GAC). The sample con¬ 
sists of 94,332 FGK main-sequence stars with well-determined radial velocities and atmo¬ 
spheric parameters. To derive the LSR, two independent analyses are applied to the data. 
Firstly, we determine the solar motion by comparing the observed velocity distribution to 
that generated with the analytic formulism of Schonrich & Binney that has been demon¬ 
strated to show excellent agreement with rigorous torus-based dynamics modelling by Bin¬ 
ney & McMillan. Secondly, we propose that cold populations of thin disc stars, selected by 
applying an orbital eccentricity cut, can be directly used to determine the LSR without the 
need of asymmetric drift corrections. Both approaches yield consistent results of solar mo¬ 
tion in the direction of Galactic rotation, Vq, that are much higher than the standard value 
adopted hitherto, derived from Stromgren’s equation. The newly deduced values of Vq are 
1-2 kms -1 smaller than the more recent estimates derived from the Geneva-Copenhagen 
Survey (GCS) sample of stars in the solar neighbourhood (within 100 pc). We attribute the 
small difference to the presence of several well-known moving groups in the GCS sample 
that, fortunately, hardly affect the LSS-GAC sample. The newly derived radial ( Uq ) and ver¬ 
tical (Wq) components of the solar motion agree well with the previous studies. In addition, 
for all components of the solar motion, the values yielded by stars of different spectral types 
in the LSS-GAC sample are consistent with each other, suggesting that the local disk is well 
relaxed and that the LSR reported in the current work is robust. Our final recommended LSR 
is, (Uq, Vq, Wq) = (7.01 ±0.20, 10.13 ± 0.12, 4.95 ± 0.09) kins” 1 . 

Key words: Galaxy: kinematics and dynamics-Galaxy: solar neighbourhood-Galaxy:stellar 
content-Galaxy: fundamental parameters. 


1 INTRODUCTION 

The local standard of rest (LSR) is defined a point in space that 
is moving on a perfectly circular orbit around the Galactic cen¬ 
tre at the solar position. An accurate estimate of the LSR, namely 
measuring the peculiar velocities of the Sun, is a starting point for 
Galactic kinematic and dynamic studies since all measurements are 
made relative to the Sun. It is needed for the derivations of basic pa¬ 
rameters such as the Oort constants, which describe the local spa¬ 
tial variations of the stellar velocity field (Feast & Whitelock 1997, 
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Oiling & Dehnen 2003), as well as the orbital eccentricity distri¬ 
butions of different stellar populations which can be used to con¬ 
strain the various scenarios of Galactic disk formation (Sales et al. 
2009, Wilson et al. 2011, Lee et al. 2011). As a fundamental task of 
Galactic astronomy, an accurate determination of the three Galactic 
space velocity components of the solar motion with respect to the 
LSR, Uq,Vq and Wq, has been a long-standing challenge. Similar 
to the ones given in Francis & Anderson (2009) and Co§kunoglu et 
al. (2011), Table 1 presents an updated summary of determinations 
(including results from the current work) of the solar motion from 
the literature of over a century. As the Table 1 shows, in spite of 
utilizing a variety of methods and tracers, the measurements have 
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hitherto failed to converge, especially for the component in the di¬ 
rection of Galactic rotation (Vq). 

In principle, the radial ( Uq ) and vertical (Wq) components 
of the solar motion can be directly obtained by taking the negative 
of the mean heliocentric velocities of various stellar populations in 
the solar-neighbourhood that are well relaxed. The situation for the 
component in the direction of Galactic rotation (Vq) is more com¬ 
plicated, due to the non-negligible values of asymmetric drift, V a , 
of stars with respect to the LSR - motions that are related to the ve¬ 
locity dispersions of the individual stellar populations of concern. 
The classical method to derive V© is based on Stromberg’s equation 
(Stromberg 1946), which assumes a linear correlation between the 
mean negative heliocentric azimuthal velocity, 14 = 14 + Vq, and 
the squared radial velocity dispersion, or, of any stellar sample. 
Then a linear extrapolation to 4 = 0 (kms x ) 2 for an idealized 
population yields 14 = Vq. Theoretically, for any well-relaxed 
stellar population, the Jeans equation, i.e. Eq. (4.228) of Binney 
& Tremaine (2008), links the asymmetric drift, V a , with the ra¬ 
dial velocity dispersion and the square bracket term (radial scale 
lengths and the properties of velocity ellipsoid) in axisymmetric 
equilibrium. The physical underpinning of Stromberg’s equation is 
that the square bracket term is a constant for the different stellar 
populations and the stellar populations used are in axisymmetric 
equilibrium. While the assumptions sound reasonable, the canon¬ 
ical value of Vq = 5.25 kms -1 , derived by Dehnen and Bin¬ 
ney (1998, hereafter DB98) based on Stromberg’s equation, has 
recently been found to be ~ 6 — 7kms -1 smaller than those de¬ 
termined using open clusters (Piskunov et al. 2006) and masers in 
massive star forming regions (Reid et al. 2009; Bobylev & Bajkova 
2010; McMillan & Binney 2010) as tracers. The possible causes 
of discrepancies have been addressed recently by a number of in¬ 
vestigators. Schonrich et al. (2010) show that Stromberg relation 
could underestimate Vq due to the presence of a metallicity gradi¬ 
ent, which breaks the linear relation between the mean rotation ve¬ 
locity and the squared velocity dispersion of stars binned by colour. 
This implies that the square bracket term may vary from one popu¬ 
lation to another population as represented binned by the individual 
colour bins. In addition, the axisymmetric equilibrium of Galactic 
potential, a basic assumption underlying Stromberg’s equation as 
mentioned above, can be broken due to the presence of significant 
non-axisymmetric structures in our Galaxy (e.g. the spiral arms, the 
central bar; Francis & Anderson 2009; Binney 2010). Apart from 
Vq, estimates of the radial and vertical components of the solar 
motion. Uq and Wq, can also be biased by the presence of non- 
axisymmetric structures. 

To minimize the effects of asymmetric drifts when estimat¬ 
ing Vq, two different methods have been developed that avoid the 
assumption of a constant square bracket term for the different stel¬ 
lar populations underlying Stromberg’s equation. One is to deter¬ 
mine the solar peculiar velocity from the offset of the modelled 
velocity distribution that best matches the observed data (Binney 
2010; Schonrich et al. 2010; Schonrich & Binney 2012, hereafter 
SB 12). Alternatively, the LSR can also be determined utilizing cold 
populations of thin disc stars as tracers, under the assumption that 
they suffer from negligible asymmetric drifts (Francis & Anderson 
2009; Co§kunoglu et al. 2011). Compared to the method based on 
Stromberg’s equation, the latter approach has the advantage that it 
is simple and straight forward, and the results do not depend on the 
assumption that underpins the former method. The key to a success¬ 
ful application of this technique is to select and define a well rep¬ 
resenting sample of cold populations of thin disc stars. It should be 
emphasised that all the approaches outlined above assume axisym- 


Figure 1 . Pseudo-color HR diagram of ~ 0.7 million stars selected from the 
LSS-GAC DR1 that have spectral signa-to-noise ratios at 4650 A, S/N(4650 
A) 7 ' 10 and robust radial velocity and stellar parameter determinations. 
The colorbar shows the stellar density on a logarithmic scale. The blue box 
delineate the parameter space used to select FGK main-sequence stars ana¬ 
lyzed in the current work. 

metric equilibrium. Thus it is essential to select a sample of stars 
that are least affected by those known non-axisymmetric structures. 
More discussions of this issue will be made below. 

In this paper, we re-estimate the peculiar velocity of the Sun 
with respect to the LSR by applying both the aforementioned meth¬ 
ods to a sample of stars selected from the Large Sky Area Multi- 
Object Fiber Spectroscopic TelescopiQ (LAMOST; also named 
the Guoshoujing Telescope) Spectroscopic Survey of the Galactic 
Anti-center (LSS-GAC, Liu et al. 2014). For the velocity distribu¬ 
tion fitting (VDF) method, we use the analytic formulism of SB 12, 
which has been shown to reproduce the results of rigorous torus- 
based dynamics modelling (Binney & McMillan 2011) to high fi¬ 
delity. As for the method based on cold populations of thin disc 
stars (CTDS), we use the orbital eccentricity, in replacement of the 
eccentricity vector adopted by Francis & Anderson (2009), as a 
natural indicator of stellar population to select cold populations of 
thin disc stars. The sample selected from the LSS-GAC comprises 
FGK main-sequence stars that satisfy the following criteria: (i) sur¬ 
face gravities 3.8 < logg < 5.0 dex; (ii) effective temperatures 
4, 200 < T e fj < 6, 800 K; (iii) distances less than 600 pc; (iv) 
total Galactic space velocity errors less than 15 km s -1 . The first 
two constraints ensure that the stellar parameters of selected stars 
are well determined with the current LAMOST stellar parameter 
pipeline (cf. Section 3). The third constraint is to minimize the ef¬ 
fects of large-scale streaming motions, induced by, for example, 
non-axisymmetric structures. The fourth constraint is to ensure ro¬ 
bust Galactic space velocities (cf. Section 4). To apply the CTDS 
method, we apply an additional cut to the sample: orbital eccen¬ 
tricities less than 0.13. The constraint is pivotal to the method (cf. 
Section 6). Stars of very low orbital eccentricities move in near¬ 
circular orbits with negligible asymmetric drifts and thus can be 
used to determine Vq without the need of correcting for their ef¬ 
fects. 

The paper is organized as following. Section 2 presents a brief 

1 LAMOST is a 4 meter quasi-meridian reflecting Schmidt telescope 
equipped with 4000 fibers, each of an angular diameter of 373 projected 
on the sky, distributed in a circular field of view (FoV) of 5° in diameter 
(Cui et al. 2012). 
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Table 1 . Measurements of the LSR in the literatures and from the current work 


Source 

Data 

Uq 

v© 

w© 



(km s 1 ) 

(km s x ) 

(km s —1 ) 

This study (2014) 

LSS-GAC DR 1 

7.01 ±0.20 

10.13±0.12 

4.95±0.09 

Bobylev & Bajkova (2014) 

Young objects 

6.00±0.50 

10.60±0.80 

6.50±0.30 

Co§kunoglu et al. (2011) 

RAVE DR3 

8.50±0.29 

13.38±0.43 

6.49±0.26 

Bobylev & Bajkova (2010) 

Masers 

5.50±2.2 

11.00±1.70 

8.50±1.20 

Breddels et al. (2010) 

RAVE DR2 

12.00±0.60 

20.40±0.50 

7.80±0.30 

Schonrich et al. (2010) 

Hipparcos 

11 in+°- 69 

11.lU_o 75 

12 24”*~ 0 ' 47 
lz - Z4 -0.47 

7 

'■ z;? -0.36 

Reid et al. (2009) 

Masers 

9.0 

20 

10 

Francis & Anderson (2009) 

Hipparcos 

7.50±1.00 

13.50±0.30 

6.80±0.10 

Bobylev & Bajkova (2007) 

F & G dwarfs 

8.70±0.50 

6.20±2.22 

7.20±0.80 

Piskunov et al. (2006) 

Open clusters 

9.44±1.14 

11.90±0.72 

7.20±0.42 

Mignard (2000) 

K0-K5 

9.88 

14.19 

7.76 

Dehnen & Binney (1998) 

Hipparcos 

10.00 ±0.36 

5.25±0.62 

7.17±0.38 

Binney et al. (1997) 

Stars near South Celestial Pole 

11.00±0.60 

5.30±1.70 

7.00±0.60 

Mihalas & Binney (1981) 

Galactic Astronomy (2nd Ed.) 

9.00 

12.00 

7.0 

Homann (1886) 

Solar neighborhood stars 

17.40±11.2 

16.90± 10.90 

3.60±2.30 
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Figure 2. Distributions of the integrated E(B — V) reddening in the solar neighbourhood of 600 pc for the current LSS-GAC DR1 sample of FGK stars in 
the Galactic X-Y (left panel) and X-Z (right panel) planes. The colorbars indicate values of E(B — V) in units of magnitude. 


description of the LSS-GAC. The data are presented in Section 3. 
Section 4 describes the Galactic space velocities of sample stars and 
uncertainty control. Section 5 and 6 presents respectively the VDF 
and CTDS methods and their results. The results are discussed in 
Section 7. Conclusions are given in the Section 8. 


2 LSS-GAC 

The LSS-GAC is a major component of the on-going LAMOST 
Galactic surveys, aiming to collect under dark and grey lunar con¬ 
ditions optical (AA3800-9000), low resolution (R ~ 1, 800) spec¬ 
tra for a statistically complete sample of over a million stars of 
all colors and of magnitudes 14.0 ^ r < 17.8 mag (18.5 mag 
for limited fields), in a continuous sky area of ~ 3, 400 sq.deg., 
centred on the Galactic anticentre (GAC), covering Galactic lon¬ 
gitudes 150 < l < 210° and latitudes |6| < 30°. Over 1.5 
million very bright stars, brighter than ~14.0mag., selected with 
a similar target selection algorithm, will also be observed, utilizing 
bright lunar conditions. The survey will deliver spectral classifica¬ 


tions, values of stellar radial velocity V r and atmospheric param¬ 
eters (effective temperature T e g, surface gravity log g, metallicity 
[Fe/H]) for about 3 million Galactic stars. The survey will sample 
main-sequence stars out to a few kpc, and tens of kpc for giants. 
Combined with distances, e.g. the spectrophotometric distances of 
Yuan et al. (2015) and the proper motions from other surveys, the 
LSS-GAC provides a huge data set that allows one to study the stel¬ 
lar populations, kinematics and chemistry of the Galactic thin/thick 
disks and their interface with the halo in multi phase space (three 
dimensional space position and velocity, and metallicity). The Pilot 
and Regular Surveys of LSS-GAC were initiated in September of 
2011 and 2012, respectively. The Regular Survey is expected to last 
for 5 years. More details about the survey, including the scientific 
motivations, target selections and data reduction, can be found in 
Liu et al. (2014) and Yuan et al. (2015). 
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Figure 3. Distribution of distances of FGK main-sequence stars in the LSS- 
GAC DR 1 within 1 kpc of the Sun. The shaded area shows stars within 
600 pc analyzed in the current work. 


3 DATA 

The data used in the current study are from the first data release 
of value-added catalogs of the LSS-GAC (hereafter DR1; Yuan 
et al. 2015). There are about 0.7 million stars in the DR1 with 
a spectral signal-to-noise ratio per pixel (~ 1.07 A) at 4650 A, 
S/N (4650 A) > 10, and with robust determinations of stellar 
radial velocity V r and atmospheric parameters ( T e g, log g and 
[Fe/H]) determined with the LAMOST Stellar Parameter Pipeline 
at Peking University (LSP3; Xiang et al. 2015a,b). The HR dia¬ 
gram of the stars is presented in Fig. 1. Parameters of FGK main- 
sequence stars are best determined with the LSP3. Comparisons 
of results deduced from LAMOST multi-epoch duplicate observa¬ 
tions, with results for common targets as given by the SDSS DR9 
(Ahn et al. 2012), RAVE (Steinmetz et al. 2006) and APOGEE 
(Ahn et al. 2013) surveys and from the PASTEL compilation 
(Soubiran et al. 2010), as well as applying the Pipeline to open 
and globular clusters, show that for FGK dwarfs the LSP3 has 
achieved an accuracy of 5.0kms _1 , 150 K, 0.25 dex, 0.15 dex 
for V T , T e ff, log g and [Fe/H], respectively. We apply a surface 
gravity cut, 3.8 < logg < 5.0 dex and an effective temperature 
cut, 4, 200 < T e ff < 6, 800 K, to select the sample stars of FGK 
dwarfs. In total, ~ 0.4 million stars are selected, as rounded by 
the blue box in Fig. 1. To test the robustness of the methods used 
in the current work, we further divide the total sample into three 
sub-samples, consisting respectively of stars of spectral type F, G 
and K, based on the effective temperatures of the stars: F-type of 
6, 000 ^ T ef f ^ 6, 800 K, G-type of 5, 300 < T e ff < 6,000 K and 
K-type of 4, 200 ^ T e ff < 5, 300 K and apply those methods to 
each sub-sample. 


3.1 Extinctions and distances 

Since most of the sample stars fall in the Galactic disc at distances 
beyond 100 pc (see Fig. 3), corrections for the interstellar extinc¬ 
tion by dust grains are essential for robust distance determinations. 
Here we adopt estimates of extinction for individual stars as de¬ 
duced from the photometric colours with the “standard pair” tech¬ 
nique (Stecher 1965; Massa et al. 1983; Yuan et al. 2013, 2015). A 
comparison with the values given by the extinction map of Schlegel 
et al. (1998) for high Galactic latitude regions shows the technique 


has achieved a precision of about 0.04 mag in E(B — V (0- Fig. 2 
shows the distributions of integrated E(B — V) out to a distance of 
600 pc in the X-Y and X-Z Galactic planes for the sample stars. 
As expected, the reddening increases with distance. The local bub¬ 
ble of very low extinction is also clearly visible in quadrants of 
Y < 0 out to 600 pc. 

For distances, geometric parallaxes are not available for the 
majority of stars targeted by the LSS-GAC. We have thus devel¬ 
oped an empirical relation between the absolute magnitudes and 
the atmospheric parameters (T e ff, log g and [Fe/H]), based on the 
MILES (Sanchez-Blazquez et al. 2006) spectral library (Yuan et 
al. 2015). The relation is based on the absolute magnitudes derived 
from the Hipparcos parallaxes (typical errors ^10 per cent) and the 
atmospheric parameters compiled and homogenized by Cenarro et 
al. (2007), mostly deduced from high spectral resolution measure¬ 
ments. We find that distances thus derived for FGK dwarfs are bet¬ 
ter than 15 per cent and show no obvious systematic errors. The 
distribution of distances of FGK dwarfs within 1 kpc in the DR1 
is presented in Fig. 3. For the purpose of current work, only stars 
within 600 pc, as shown by the shaded area are analyzed. 

Fig. 4 shows the spatial distributions of the selected stars in 
the Galactic planes for a right-handed Cartesian coordinate system 
centred on the Galactic centre, where the X-axis passes through 
the Sun and points towards the Galactic centre, the Y-axis is in 
the direction of Galactic rotation and the Z-axis points towards the 
North Galactic Pole (NGP). The Sun is located at (X, Y, Z ) = 
(—Ro, 0, 0), where Ro = 8 kpc. Most of the stars fall in the anti¬ 
centre direction, with the rest coming from the Very Bright (VB) 
plates observed under bright lunar conditions. 

3.2 Radial velocities and proper motions 

As mentioned above, for FGK dwarfs radial velocities yielded by 
the LSP3 are accurate to ~ 5 km s -1 . However, we find that for 
common stars targeted by both surveys, the LSS-GAC velocities 
yielded by the LSP3 have an offset of 3.1 km s -1 relative to those 
of APOGEE, the latter have a measurement accuracy better than 
0.15 km s" 1 . Hence, we have applied an offset correction to all 
radial velocities yielded by the LSP3 (cf. Xiang et al. 2015b for a 
detailed description). 

The proper motions of our sample stars are taken from the 
Fourth United States Naval Observatory (USNO) CCD Astrograph 
Catalog (UCAC4; Zacharias et al. 2013), an all-sky astrometric cat¬ 
alog archiving over 113 million objects. Among them, over 105 
million stars have proper motions. The catalog is completed to 
R ~ 16 mag. Typical random and systematic errors of the UCAC4 
proper motions are 4 and 1- 4 mas yr -1 , respectively (Zacharias et 
al. 2013). 

For our sample, a systematic error of a few mas yr -1 in proper 
motions could induce a bias of several km s -1 in the determina¬ 
tion of LSR, especially for the Vq and W q components. To cor¬ 
rect for the systematic errors of UCAC4 proper motions, we have 
cross-matched the UCAC4 with the 2nd release of Large Quasar 
Astrometric Catalog (LQAC, Souchay et al. 2012). After exclud¬ 
ing quasars with poor proper motion determinations, this yields a 
total of 1,700 quasars distributed over the whole sky. Corrections 

2 For targets without optical photometry, the estimated uncertainties of ex¬ 
tinction increase to about 0.1 mag in E(B — V) (Yuan et al. 2015). 

3 Quasars, being distant extragalactic point sources, have essentially zero 
proper motions and are thus the best objects to estimate and correct for 
uncertainties of stellar proper motion measurements. 
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Figure 4. Spatial distributions of the sample stars in the X - Y (left panel) and X — Z (right panel) planes. The logarithmic stellar number densities in units 
of number per kpc 2 are indicated by the side colorbars. 


Table 2. Coefficients of a two-dimensional second-order polynomials fit to the systematic errors of proper motions of UCAC4 as a function of the Right 
Ascension a and Declination 8. 


Components 

a. 0 

ai 

«2 

a 3 

CI 4 

«5 

/io: cos S 

1.29 

-0.0137 

3.27 x 10“ 5 

-4.90 x 10“ 5 

-0.0106 

5.01 x 10" 5 

V5 

-1.82 

3.68 x 10“ 3 

-2.91 x 10“ 6 

-5.57 x 10“ 5 

-4.70 x 10“ 3 

2.35 x 10" 4 


Note: The fitting function has the form: /r = ao + aia + a2Ct 2 + 7/307) + 0.4 /) + <15 < 5 2 , where //.©, in mas yr —1 . 
represent the the two components of quasar proper motions, fi a cos <5 and jig, and a and <5, both in units of degree, 
are the Right Ascension and Declination. 


to account for the systematic errors of UCAC4 proper motions are 
then derived by fitting the proper motions of quasars (by default, 
quasars should have zero proper motions) with a two-dimensional 
second-order polynomials as a function of the Right Ascension and 
Declinatior0. Table 2 presents the fit coefficients. No magnitude 
and color dependences of the systematic errors of proper motions 
of UCAC4 are found for quasars of R magnitudes between 12 and 
16 mag and B — V color between 0.1 and 1.4 mag. UCAC4 proper 
motions corrected for the systematic errors, as calculated from the 
above fits, are used for all calculations given below. 


4 GALACTIC SPACE VELOCITY AND UNCERTAINTY 
CONTROL 

From the distances calculated with the empirical relation described 
in Section 3., the LSP3 radial velocities and the UCAC4 proper 
motions, with the latter two corrected for the systematic errors as 
discussed in Section 3.2, we calculate the total Galactic space ve¬ 
locities and their three components, ( U , V, W ) of the sample stars 
using the standard transformation matrices derived by Johnson & 
Soderblom (1987). U, V and W refer to respectively the velocity 
components in the X, Y and Z directions in a right-handed Carte¬ 
sian coordinate system as defined in Section 3.1. 

Uncertainties of the three velocity components, 5u, Sv and 

4 Similar corrections for the PPMXL (Roeser et al. 2010) proper motions 
have been obtained by Carlin et al. (2013). 


5w, are also calculated using the standard algorithm given by John¬ 
son & Soderblom (1987) by propagating the uncertainties of the ra¬ 
dial velocities, proper motions and distances. The uncertainties of 
total Galactic space velocities are then calculated using, 

= Su + $V + ( 1 ) 

The distributions of uncertainties of the total Galactic space 
velocities and their three components are presented in Fig. 5. The 
mean and standard deviation of uncertainties of the total Galactic 
space velocities are ( 6s ) = 10.91 km s _1 and s.d. = 8.25 km 
s -1 , respectively. Stars with 5s > 15 km s -1 are discarded, to¬ 
taling 30,829 stars, about 25% of the sample. After applying the 
cut, the average uncertainties of the three velocity components for 
the remaining stars of the sample are reduced to about 5 km s _ , 
with a standard deviation of ~ 2 km s -1 . The velocities of the re¬ 
maining sample stars are thus sufficient accurate to make a robust 
determination of the LSR. 

5 DETERMINING V© BY VDF 
5.1 Bulk motion effects 

As pointed out earlier, determining the solar peculiar velocity in 
the Galactic rotation direction, V©, is not as straightforward as for 
those in the other two directions (( 7 © and IV©), due to the presence 
of significant asymmetric drift in the azimuthal direction, V a . The 
classical linear Stromberg’s equation was proposed to solve this 
problem (e.g. DB98). However, as pointed out by Schonrich et al. 
(2010), Stromberg’s equation can be broken due to the presence of 
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Figure 5. Histograms of Galactic space velocity uncertainties. The top 
panel shows the total Galactic space velocities and the rest the three com¬ 
ponents. Stars with a total Galactic space velocity error 5s > 15 km s — 1 
are removed from the sample. The histogram distributions of the remaining 
stars are represented by the shaded areas. The mean and standard deviation 
of errors of the total velocities marked in the top panel are calculated be¬ 
fore applying the cut, whereas those in the rest panels for the three velocity 
components are calculated after applying the cut. 

a metallicity gradient that breaks the linear relation between asym¬ 
metric drift and radial velocity dispersion of stars binned by colour. 
To overcome the obstacle, a technique has been developed to de¬ 
rive Vq from the offset of the velocity distribution of a dynamical 
model that best matches the measurements. Using this technique, 
Binney (2010) constructs an azimuthal velocity distribution based 
on actions estimated from the adiabatic invariance of the local sam¬ 
ple of stars of the Geneva-Copenhagen Survey (GCS, Nordstrom et 
al. 2004; Holmberg, Nordstrom & Andersen 2009), and revises Vq 
upwards to 11 km s -1 . Similarly, by chemodynamically modeling 
the GCS sample, Schonrich et al. (2010) conclude that DB98 un¬ 
derestimated Vq by ~ 7 km s _1 . As described earlier, the VDF 
method is also founded on the assumption of axisymmetric equilib¬ 
rium. On the other hand, it has been known that non-axisymmetric 
structures (e.g. the spiral arms, the central bar) could distort the 
observed velocity distribution (Dehnen 1999; De Simone, Wu & 
Tremaine 2004; Antoja et al. 2009; Minchev et al. 2009). For the 
moment, the only way to minimize the dynamical effects of those 
non-axisymmetric structures is to carefully select a sample of stars 
least affected by those structures. Unfortunately, a number of mov¬ 
ing groups have been found in the ( U , V) velocity space in the GCS 


sample (e.g. Dehnen, 1998), which is restricted to the immediate 
vicinity of the Sun (within 100 pc) and consists of just over ten 
thousand stars of a variety of spectral types. As shown in Fig. 6, 
the distribution of the V velocity component of GCS stars is seri¬ 
ously distorted by four well-known moving groups, especially the 
Pleiades and/or Hyades streams. Unless the effects of bulk motions 
of those streams could be properly modeled, the aforementioned 
determinations of Vq based on the GCS sample could be signifi¬ 
cantly biased. Unfortunately, the origin of those moving groups is 
still on debate (e.g. Antoja et al. 2011), making it difficult to either 
incorporate them in the dynamical model or remove their effects 
by subtracting the overdensities produced by them in the observed 
velocity distribution. Another large sample of local stars, alterna¬ 
tive to the GCS one and less plagued by bulk motions of moving 
groups, is needed to apply this powerful technique developed by 
Biney (2010) in order to accurate determine Vq without biases. 
The LSS-GAC sample of local stars presented in the current work 
seems to exactly fill the bill. As shown in Fig. 7, the distribution of 
Vtf, velocity components of our LSS-GAC local sample of stars is 
on the whole quite smooth. Clearly, compared to the GCS sample, 
the LSS-GAC sample presented here is much less affected by the 
bulk motions, and thus should enable us to obtain a more robust 
estimate of Vq than possible with the GSC sample. 

5.2 Determining Vq by fitting the velocity distribution of the 
LSS-GAC sample 

To determine Vq, the model azimuthal velocity distribution is com¬ 
pared to the observed one of the LSS-GAC sample. The model 
azimuthal velocity distribution is generated using the analytic for¬ 
mula derived by SB12. As argued by SB 12, this analytic formula 
provides excellent fit to the azimuthal velocity distribution yielded 
by rigorous torus-based dynamics modelling (Binney & McMillan 
2011) and reproduces the observed distribution of the GCS local 
sample of stars, which have accurate space velocity measurements. 
We refer the reader to SB 12 for a detailed discussion about this for¬ 
mula. From the formula, the distribution of Vcf, at a given Galactic 
position ( R , z) is given by, 

n{V*\R,z) =Ue-^- R ^ /R -‘^£. 

a ( 2 ) 
x exp(- A$ad ) Z(z, i?g - R ), 

where N is a normalizing factor. R g is the guiding-center radius 
given by, R g = RV^/Vc (where V c is the circular velocity as a 
function of R g , namely, the rotation curve. In this work we adopt 
a flat rotation curve of constant value 220 km s -1 ). Rq = 8 kpc is 
the Galactocentric distance of the Sun, Rd the disc scale length, K 
a factor that can be numerically determined [cf. Eq (12) of SB12]. 
The term cr, being a function of R g , is given by 

a(R g ) = a 0 e~ (R ~ Rs)/R ' r , (3) 

where op is the local velocity dispersion, R a the scale length on 
which the velocity dispersion varies. AT> a d is the so-called “adia¬ 
batic potential” that includes the effects of radial motion induced 
by adiabatic invariance of vertical motion [cf. Eq (33) of SB 12 for 
detail], f(z,R g — R) is the z factor (cf. Section 3.1 in SB 12), which 
introduces other two parameters, and the disc scale length h z . 

5 SB 12 assume that the vertical force K z is proportional to z a ~ 1 . Most 
of the LSS-GAC sample stars are near the Galactic midplane and we adopt 
a = 1.5 in the current work. 
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Figure 6. Distribution of the V velocity components of the GCS sample of 
local stars. Four well-known moving groups are clearly visible. 


Table 3. Parameters adopted for the two components employed to fit the 
velocity distribution. 


Components 

R d (kpc) 

h z (kpc) 

Rv (kpc) 

Cool 

2.4 

0.2 

7.5 

Hot 

2.5 

1.0 

15.0 


We fit the V,f, distribution of our LSS-GAC sample by two 
components, a cool component representing the thin disc and an¬ 
other hot component representing the thick disc. The distribution 
of each of the two components are described by the above formula. 
The parameters adopted for the two components are listed in Table 
3. To fit the data, each component is allowed to have two free pa¬ 
rameters: the normalizing constant M and the local velocity disper¬ 
sion <To. The combined distribution of V^ of the two components is 
then calculated for a typical position (R = 8.31 kpc, z = 0.04 kpc) 
of our sample stars. The distribution is then shifted by an offset, 
which gives the value of Vq, in order to fit the observed distribu¬ 
tion of the LSS-GAC sample. Fig. 7 shows the best fit performed for 
the velocity range, 130 < < 275 km s^ 1 , on a linear and loga¬ 

rithmic scale. The two component fit reproduces the observed dis¬ 
tribution remarkably well. For the whole sample, the best fit yields 
Vq = 9.75 ± 0.19kms _1 . The local velocity dispersions of the 
cool and hot components are 21.61 and 35.47 km s -1 , respectively. 
The analysis is repeated for the three sub-samples, consisting re¬ 
spectively of stars of spectral type F, G and K, selected based on 
the values of effective temperature as described in Section 3. 

The typical position (R, z), the number of stars N, as well 
as the value of Vq deduced from the individual sub-samples, are 
listed in Table 4, along those for the whole sample. Table 4 shows 
that values of Vq yielded by the individual sub-samples are con¬ 
sistent with each other suggesting the robustness of the results. The 
analyses also show that the local velocity dispersions of F-type stars 
are smaller than those of G- and K-type. The result can be easily 
explained by the fact that stars of F-type are younger on average 
compared to those of G- and K-type, thus have smaller velocity 
dispersions. The same effect also explains why the value of Vq de¬ 
termined from F-type are slightly smaller than that deduced from 
G- or K-type stars. 


6 DETERMINING THE SOLAR MOTION BY CTDS 

6.1 Orbital eccentricities 

In addition to the technique of VDF, another method to determine 
the LSR is using only stars from the cold populations, thin disc that 
suffer from negligible asymmetric drifts. The method is simple and 
straightforward. Co§kunoglu et al. (2011) analyze the populations 
of a RAVE sample from which they select a sub-sample of thin 
disc stars based on their Galactic space velocity distribution, as 
proposed by Bensby et al. (2003, 2005), to derive the LSR. The re¬ 
sult again confirms that previous estimates of Vq from Stromberg’s 
equation are grossly underestimated by 6-7 km s -1 . However, the 
Galactic space velocity distribution is not an efficient way to dis¬ 
entangle thin and thick disc stars, as shown by Binney (2010). As 
a result, the sample of thin disc stars selected by Co§kunoglu et al. 
(2011) may be biased. Here we use the orbital eccentricity [defined 
by Eq. (2)] as a natural indicator to define and select cold popula¬ 
tions of thin disc stars. Stars of low orbital eccentricities move in 
nearly circular orbits with negligible asymmetric drifts, whereas the 
orbits of those of stars with large eccentricities deviate from circu¬ 
lar ones and thus possess significant asymmetric drifts. In addition, 
by applying a stringent eccentricity cut, stars with significant bulk 
motions can also be screened out (Francis & Anderson 2009). By 
applying an eccentricity cut to the current sample, we expect to de¬ 
termine the LSR straightforwardly and robustly. For this purpose, 
we use an eccentricity cut of 0.13, the peak value of thin disc stars 
(Lee et al. 2011). As shown by Bobylev and Bajkova (2014), the 
orbital eccentricities of young O to B2.5 stars, used by them to de¬ 
rive the LSR range from 0 to 0.13. The mean orbital eccentricity of 
Galactic open clusters is also found to be around 0.13 (Piskunov et 
al. 2006). All those findings suggest that 0.13 is a reasonable eccen¬ 
tricity cut for the purpose of selecting thin disc stars. Furthermore, 
we will show in Section 7 that the asymmetric drifts of cold popu¬ 
lations selected using this orbital eccentricity cut are negligible. 

The eccentricity of a star is defined by, 

77apo 77peri ... 

6 = 7?— tj r~ ’ (4) 

where i? apo (-R per i) is the maximum (minimum) Galactocentric 
distance reached by the star in its orbit. To calculate the eccen¬ 
tricities of stars of our FGK main-sequence sample, we integrate 
their orbits assuming a Galactic potential given by Gardner & Flynn 
(2010). In their model, the Milky Way consists of a Miyamoto & 
Nagai (1975) disc, a Plummer (1911) bulge and inner core, and 
a spherical logarithmic halo. The characteristic parameters can be 
found in their Table 1. 

6.2 Impacts of the initial LSR 

To calculate orbital eccentricities, one however needs to assume a 
LSR. In our approach, we first calculate orbital eccentricities with 
an initial guess of LSR and then iterate the entire LSR determina¬ 
tion process until all components of the solar motion converge to a 
degree less than the fitting errors (cf. Section 6.3). In this approach, 
it is thus of utter importance to check and ensure that the final val¬ 
ues of LSR thus derived are not sensitive to the initial values of 
the LSR assumed. We have thus carried out extensive test of the 
impacts of initial LSR on the final results. As shown in Table 5, 
we have used three vastly different sets of estimates of the LSR 
from the literature as our initial guess. As pointed out earlier, the 
Uq and Wq components of the solar motion are hardly affected by 
the asymmetric drifts and straightforward to estimate. It is thus not 
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Figure 7. The V' 0 velocity distribution of the LSS-GAC sample of FGK main-sequence stars are fitted by the combination of two components: a cool 
component for the thin disc and a hot component for the thick disc. Each component, characterized by distribution function Eq. (2), has two free parameters: 
a normalization constant A f and the local velocity dispersion erg. To fit the data, the model distribution is shifted by a constant offset, which represents the 
peculiar velocity of the Sun, i.e. Vq. The left panel shows the model distributions and data points on a the linear scale, while in the right panel show the same 
distributions and data on a logarithmic scale, which reveals better the wings of the distributions. 


Table 4. The solar motion in the Galactic rotation direction with respect to the LSR Vq determined by VDF. 


Sample 

T e « 

(K) 

v 0 

(km s *) 

R 

(kpc) 

Z 

(kpc) 

■"* on 
tz> O 

1 £ 

hot 

°0 

(km s 1 ) 

N 

Whole 

4200 T eff ^ 6800 

9.75 ± 0.19 

8.31 

0.03 

21.61 

35.47 

94332 

F type 

6000 «: T eS ^ 6800 

9.00 ±0.29 

8.38 

0.04 

19.02 

26.85 

24220 

G type 

5300 ^ T ef f < 6000 

10.03 ± 0.21 

8.33 

0.05 

23.21 

38.90 

39100 

K type 

4200 sC T eff < 5300 

10.20 ±0.23 

8.24 

0.01 

22.27 

38.29 

31012 


surprising that their final derived values for the three completely 
different sets of initial guess do not differ much. It is remarkable 
that, for the Vq component which is known to suffer from signif¬ 
icant asymmetric drift effects, the final estimates that result from 
the three completely different sets of initial guess that span a wide 
range of value from ~ 5 to 20 km s _1 all converge within one per 
cent! In determining the LSR, we first calculate the orbital eccen¬ 
tricities of sample stars for each assumed set of initial values of 
LSR and then apply orbital eccentricity cut to select a sample of 
cold populations. The three components of the solar motions with 
respect to the LSR are then determined by the method described in 
the next subsection using the selected sub-sample of cold popula¬ 
tions. As shown in Table 5, the final values of LSR converge after 
4-5 iterations (the numbers of which are listed in the last column 
of the Table 5). The fact that the final values of LSR are insensitive 
to the initial guess can be understood physically. Stars belonging to 
the thick disc or halo population generally lag behind the rotation of 
thin disc stars by more than 45 km s _1 , a lag that is more than two 
times the largest value of the assumed initial ones of Vq . This im¬ 
plies that the orbital eccentricities of most thick disc or halo stars, 
calculated for the assumed initial LSR, remain too large to pass the 
orbital eccentricity cut, even for the case of initial guess of the 
largest value of Vq . Meanwhile, most stars of our sample are near 
the Galactic midplane. The sample is thus dominated by cold popu¬ 
lations of thin disc stars. Changes in the calculated eccentricities as 
induced by the different sets of initial LSR for most of those stars 
are quite small, and thus unlikely to significantly affect our selec¬ 
tion of cold populations of thin disc stars using an eccentricity cut. 
At the beginning of the iteration process, there might be some con¬ 


taminations from those hot thick disc or halo stars. However, as the 
iteration progresses and the resultant LSRs converge to true values, 
the amount of contaminations decreases. It is thus quite natural that 
after several iterations, all results converge, even for quite different 
initial guess of LSR. 

The final distribution of orbital eccentricities of our FGK 
main-sequence sample is shown in Fig. 8 (calculated as an example 
with the final LSR derived assuming the initial LSR of DB98). The 
sample consists of 60,595 thin disc stars selected by applying an 
eccentricity cut of less than 0.13. The sample is used to derive the 
LSR in next subsection. 

6.3 Corrections for differential Galactic rotation and 
determinations of the solar motions 

Considering differential Galactic rotation, the Galactic space ve¬ 
locity components of a star measured by an observer at the Solar 
position can be expressed as follows: 

U = sin a — Uq + trj , (5) 

V = (V c - C c s cos a) - Vq + tv , (6) 

W = —Wq + tw , (7) 

where a is the angle between the star and the Sun with respect to 
the Galactic centre, V c is the circular speed of the Sun, 14 s repre¬ 
sents the mean azimuthal speed of the selected stellar population 
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Table 5. Impacts of the initial LSR. 


Sources 

Uq 

(km s 1 ) 

Initial 

Vq 

(km s - 1 ) 

w 0 

(km s -1 ) 

Uq 

(km s J ) 

Final 

Vfe 

(km s 1 ) 

Wq 

(km s -1 ) 

N 

Dehen & Binney (1998) 

10.00 

5.25 

7.17 

6.94 

10.35 

4.94 

4 

Reid et al. (2009) 

9.00 

20.00 

10.00 

6.98 

10.38 

4.95 

5 

Co§kunoglu et al. (2011) 

8.50 

13.38 

6.49 

6.99 

10.39 

4.95 

4 


Table 6. Solar motion with respect to the LSR determined by the method of CTDS. 


Sample 


T e fr 


Uq 

v© 

IT© 

N 



(K) 


(km s x ) 

(km s 1 ) 

(kms- 1 ) 


Whole 

4200 

< Tefj ^ 

6800 

6.94 ±0.25 

10.35 ±0.15 

4.94 ± 0.09 

60595 

F type 

6000 

^ T e fj ^ 

6800 

7.11 ±0.32 

10.06 ±0.17 

5.22 ± 0.10 

17206 

G type 

5300 

C T eff < 

6000 

6.56 ±0.31 

10.27 ±0.18 

4.83 ± 0.12 

23722 

K type 

4200 

C T eS < 

5300 

7.27 ±0.31 

10.66 ±0.18 

4.78 ±0.11 

19667 



0.0 0.2 0.4 0.6 0.8 1.0 

Eccentricity 


Figure 8. Distribution of orbital eccentricities of sample. A total of sam¬ 
ple 60,595 stars have eccentricities less than 0.13 and they comprise the 
subsample of thin disc stars used to calculate the LSR. 



Figure 9. (/-component velocities are plotted against sin a for the LSS- 
GAC thin disc sample of 60,595 FGK dwarfs. The blue dots represent the 
mean velocities of the individual bins of sin a and the red line is a linear fit 
to the dots. 


of concern here, and tu,ty, £w are respectively the noise values in 
the three velocity component U, V and W, resulting from the com¬ 
bined effects of the velocity measurement uncertainties 5u, Sv, <5w 
and the (intrinsic) velocity dispersions au, cry, ow ■ 

In order to determine the LSR, two approaches are generally 
used to account for the effects of differential Galactic rotation in 
Eqs. (5) and (6). The first approach is to fit the observed values of 
U and Vasa function of a (Schonrich et al. 2012). Alternatively, 
one can take a flat rotation curve with a constant V c = 220 km s~ 1 , 
and then subtract the terms of differential Galactic rotation from 
the observed values of U and V in Eqs. (5) and (6). Here we should 
point out that the 600 pc distance cut discussed above is needed for 
both the approaches to work, since the cut effectively reduces the 
effects of large-scale streaming motions, such as the radial veloc¬ 
ity gradient and pattern in an extended disc (Siebert et al. 2011, 
Williams et al. 2013), induced by, for example, non-axisymmetric 
structures (e.g. central bar, spiral arms). The large number of stars 
peaking at 600 pc in our sample (Fig. 3) also ensures that a suffi¬ 
cient number of stars is available for robust LSR determination. 

The magnitude of streaming motions in U, represented by the 
term V c a sin a in Eq. (5), is of the same order of magnitude of the 
velocity noise, on the level of dozens of kms -1 for stars of the 
current sample. However, the range of variations of the similar term 
in V, i.e. 14 s cos a, is less than 1 km s _1 , and thus will be lost in 
the velocity noise. Thus the first approach can only be used to 
determine [/©. Fig. 9 plots the velocity component U against sin a 
for 60,595 thin disc stars, selected from our LSS-GAC sample of 
FGK dwarfs after imposing an eccentricity cut of 0.13. The DB98 
initial LSR is assumed. A linear fit to the data yields ( Uq , V c s ) = 
(7.03 ± 0.35, 254.00 ± 10.54) km s" 1 . The value of V c s is not 
constrained by the current sample due to the relatively small range 
of sin a covered by the sample. 

For the second approach, after correcting for differential 
Galactic rotation, only the solar peculiar velocities and the veloc¬ 
ity noise remain. The noise value of each velocity component is 
assumed to obey a Gaussian distribution (Bensby et al. 2003, 2005) 

of dispersions (,u,v,w = \J^uvw + a u vw- Thus the space 
velocities after corrections should also follow Gaussian distribu¬ 
tions, U' = U - Vc sin a ~G(-f/ 0 , Cn), V' = V - (V c - 
V c cos a) ~ G(—V 0) Cv) and W' = W ~G(-W 0 , C w)- Dis¬ 
tributions of the three velocity components, after correcting for dif¬ 
ferential Galactic rotation, for the current sample of 60,595 thin 
disc stars are presented in Fig. 10, for the case of the DB98 ini- 
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Figure 10. Distributions of the three velocity components of a sample of 60,595 cold thin disk stars, after correcting for differential Galactic rotation in radial 
and azimuthal direction (U' and V'). Also overplotted are Gaussian fits to the distributions. Note that the distribution of U' is truncated due to the lack of stars 
of high eccentricities in the current sample as a result of the orbital eccentricity cut. Estimate of the mean of the distribution, which gives —Uq, is however 
not affected by the truncation. 


Table 7. Final adopted solar motions with respect to the LSR. 


Methods 

Uq 

(km s b ) 

^0 

(km s b ) 

Wo 

(kms” 1 ) 

VDF 

- 

9.75 ±0.19 

- 

CTDS a ’* 

6.97 ±0.25 

10.37 ±0.15 

4.95 ±0.09 

CTDS b ’* 

7.08 ±0.34 

- 

- 

Adopted 

7.01 ±0.20 

10.13 ±0.12 

4.95 ± 0.09 


a After corrected differential Galactic rotation assuming a flat 
rotation curve of constant V c = 220 km s _ 1 . 
b After corrected differential Galactic rotation by fitting 
Eq. (5) (see also Fig. 9). 

The value of each component is the mean of results deduced 
from the three sets of initial LSR listed in Table 5. 

tial LSR. Also overplotted in the Figure are Gaussian fits to the 
distributions. The corrected velocity distributions, U' and V' , are 
well described by a Gaussian function. The velocity distribution 
in radial component is however truncated as a result of the orbital 
eccentricity cut applied to the sample. The estimate of the mean 
of the distribution, which yields —Uq, is however not affected by 
the truncation. The fits presented in Fig. 10 yield (Uq, Vq, Wq) = 
(6.94 ± 0.25, 10.35 ± 0.15, 4.94 ± 0.09). The value of Uq de¬ 
rived from the Gaussian fit is consistent with that yielded by the 
first approach, i.e. by fitting differential Galactic rotation, as given 
above. 


7 DISCUSSION 

The value of Vq obtained in the current work by the VDF method 
is 1-2 km s' 1 smaller than those of Binney (2010) and Schonrich 
et al. (2010). The small differences might be caused by bulk mo¬ 
tions that affect the samples used in these work, as discussed ear¬ 
lier. As Fig. 6 shows, the presence of the Pleiades and/or Hyades 
streams moves the peak of V distribution of the GCS local sample 
toward negative values, leading to their estimates of Vq being 1- 
2 km s' 1 larger than obtained with the LSS-GAC sample, which 


shows a much smoother azimuthal velocity distribution and is less 
affected by the bulk motions of stellar streams. 

For the CTDS method, we have determined the LSR based 
on two assumptions. Firstly, we assume that cold populations of 
thin disc stars selected by applying an orbital eccentricity cut suf¬ 
fer from negligible asymmetric drift effects. The asymmetric drifts 
of cold populations of thin disc stars selected from the orbital ec¬ 
centricity cut can be estimated approximately using Eq. (4.228) of 
Binney & Tremaine (2008) from the velocity dispersions (<ju,v,w), 
velocity ellipsoid, stellar density (v) and circular speed (V c ). As¬ 
sume that the tilt angle is toward the Galactic center (Binney et al. 
2014; Budenbender et al. 2014) and both u and au are exponential 
functions of R, 


v oc exp(— 



( 8 ) 


cru oc exp( — ), (9) 

-iter 

where Rd and R a are the scale lengths of v and au, respec¬ 
tively. Under those assumptions, Eq. (4.228) of Binney & Tremaine 
(2008) becomes, 



&v ± ayy ± Ra u (——|- 


2 

Ra 



( 10 ) 


where the intrinsic velocity dispersion^ of the selected thin disc 
stars are au = 18.98, av = 11.05 and aw = 9.61 kms- 1 , 
respectively. Adopting Ra = 13.70 kpc from Sharmaet al. (2014), 
Rd = 2.5 kpc, Vc = 220 kins -1 and a typical position R = 8.3 
kpc of the selected stars, we have V a = 2.56 kms -1 . If using 
instead a larger scale length Rd = 3.7 kpc from the recent work 
of Chang et al. (2011), the estimated value of V a decrease to 1.68 
kins -1 . 


6 The intrinsic velocity dispersions are estimated by subtracting the space 
velocity errors using the sample of 60,595 thin disc stars described in Sec¬ 
tion 6.2. 
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Figure 11. Same as Fig. 10 but for sub-samples of different spectral types. The numbers of stars of the individual spectral types are marked in the left panels 
of the Figure. 


The results show that our assumption that cold populations of 
thin disk stars selected with an orbital eccentricity cut of e < 0.13 
have small asymmetric drifts is a reasonable one. This is corrob¬ 
orated by the close agreement of values of Vq deduced from the 
two independent analyses, i.e. those based on the VDF and CTDS 
methods, respectively. The effects of asymmetric drifts of the cur¬ 
rent sample of thin disc stars on Vq estimation, if any, are likely to 
be quite smalfl As a further test of the assumption, we note that 


Considering the fact that the value of Vq deduced from the CTDS 


stars of different spectral types (populations) suffer from different 
levels of the asymmetric drift effects. For sub-samples of dwarfs of 
F, G and K spectral type, the level of those effects can vary from 
a few to tens of km s -1 (DB98, Schonrich et al. 2010). The fact 
that all the three of F, G and K sub-samples give consistent results 
suggest again the assumption is a reasonable one. 


method is only 0.62 km s —1 larger than that deduced from the VDF method 
(which fully accounts for the effects of asymmetric drifts), the value of V a 
of the current sample of thin disc stars is likely to be smaller than 1 km s~ 1 . 


7 
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As described earlier, the second assumption underlying our 
analysis with the CTDS method is that the sample suffers from 
little systematic kinematic effects, induced by, for example, non- 
axisymmetric structures, accretion events or asymmetric mass in¬ 
fall. As already pointed out in Section 5.1, the current sample does 
not seem to be affected significantly by any moving group. The 
validity of this assumption also seems to be supported by the con¬ 
sistent results yielded by the individual sub-samples of different 
spectral type, considering that those sub-samples probe different 
depths (because of their different intrinsic, absolute magnitudes^ 
Fig. 11 plots the results from the three sub-samples for the case of 
the DB98 initial LSR. The values of solar motion derived from the 
sub-samples are presented in Table 6. Fig. 11 and Table 6 show that 
for all the three velocity components of LSR, the values yielded 
by all the three sub-samples agree within the fitting errors. As a 
by-product, the results also suggest that the local disk seems to be 
radially and vertically stable and circularly symmetric. 

We note all the errors quoted for the estimates of LSR pre¬ 
sented here are purely formal. There are potential systematic errors 
for both the methods of VDF and CTDS. For the VDF method, po¬ 
tential biases may come from the underlying assumptions and ap¬ 
proximations, as well as from the parameters adopted to construct 
the analytical azimuthal velocity distribution. As for the CTDS 
method, we believe that the effects of asymmetric drift are minor, 
although effects on the level of ~ 1-2 km s^ 1 could be present as 
estimated above, considering that the thin disc stars selected do not 
have exactly zero eccentricity, i.e. they do not move in perfect cir¬ 
cular orbits. Although there might be a variety of potential sources 
of systematic errors affecting the results from both methods, it is 
reassuring that both methods yield consistent values of Vq within 
2a. Our final, adopted values of LSR are presented in Table 7, ob¬ 
tained by taking the mean of individual estimates, weighted by the 
inverse square of the formal errors of individual estimates. 


8 CONCLUSIONS 

The traditional method based on Stromberg’s equation fails to pro¬ 
duce accurate estimates of the LSR, especially for the Vq com¬ 
ponent, due to the presence of kinematic and/or metallicity biases 
addressed recently by a number of investigators. In this work, two 
alternative and independent methods have been used to derive the 
LSR using a completely new sample of stars selected from the LSS- 
GACDR1. 

The first method derives V© by fiting the observed azimuthal 
velocity distribution with an analytical distribution function which 
has been shown to reproduce well reults from rigorous torus-based 
dynamics modeling. We obtain Vq = 9.75 ± 0.19kms _1 . The 
value is about 1-2 kms' 1 smaller than those derived by Binney 
(2010) and Schonrich et al. (2010) by applying the same method 
to the GCS local sample. We suggest that the small differences are 
caused by the presence of stellar streams that significantly plagues 
the GCS sample, but to a much less extent the LSS-GAC DR1 sam¬ 
ple. 

The second method derives the LSR directly using cold popu¬ 
lations of thin disc stars as tracers. The method relies on the defini¬ 
tion of a proper sample of cold populations of thin disc stars. In this 

8 Here, we are not arguing that cold populations are less affected by those 
non-axisymmetric structures than any warm tracer. We just point out, by 
good fortune or not, that the current sample of cold populations seems to 
happen to be hardly affected by those structures. 


paper, we propose that the orbital eccentricity serves as a natural 
quantity to select cold populations of thin disc stars. After applying 
an orbital eccentricity cut, two approaches are adopted to correct 
for the effects of differential Galactic rotation and determine the 
LSR. The first approach is applicable to the Uq component only. 
The analysis yields Uq = 7.03 ± 0.35 km s -1 . The second ap¬ 
proach, under the assumption of a flat rotation curve and negligible 
asymmetric drifts, is based on Gaussian fitting to the measured ve¬ 
locity distributions, after correcting for the effects of differential 
Galactic rotation. Values of Uq derived from both approaches are 
consistent with each other and in agreement with that deduced from 
maser measurements (Bobylev & Bajkova 2010, 2014). From the 
second approach, we find Vq = 10.35 ± 0.15 km s -1 , which is 
consistent within 2a with the value deduced by the method of VDF. 
Our newly derived value of Wq = 4.94 ± 0.09 km s -1 is very 
close to the recent estimate of Schonrich (2012). Both for VDF 
and CTDS methods, the whole sample is further divided into sub¬ 
samples of F, G and K spectral types. Results derived from those 
sub-samples all yield consistent results, as listed in Table 4 and 6. 
The results corroborate the robustness of the LSR determined by 
both methods. 

Benefited from the LSS-GAC DR1 and the UCAC4 proper 
motion catalog, we have built the largest local sample of FGK 
main-sequence stars within 600 pc of the Sun with accurate dis¬ 
tances and velocities. By controlling the random errors discussed 
in Section 4 and correcting systematic errors both in radial veloci¬ 
ties and in proper motions as discussed in Section 3.2, the sample 
has enabled us to determine the LSR with the smallest uncertainties 
hitherto (see Table 1). Our final recommended values of LSR are: 
(Uq, Uq, Wq) = (7.01 ± 0.20, 10.13± 0.12, 4.95 ± 0.09) km s” 1 
(see Table 7). 
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